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Abstract 

Fermionic linear optics is efficiently classically simulatable. Here it is shown that the set 
of states achievable with fermionic linear optics and particle measurements is the closure of a 
low dimensional Lie group. The weakness of fermionic linear optics and measurements can 
therefore be explained and contrasted with the strength of bosonic linear optics with particle 
measurements. An analysis of fermionic linear optics is used to show that the two-qubit match- 
gates and the simulatable matchcircuits introduced by Valiant generate a monoid of extended 
fermionic linear optics operators. A useful interpretation of efficient classical simulations such 
as this one is as a simulation of a model of non-deterministic quantum computation. Problem 
areas for future investigations are suggested. 

1 Introduction 

It is conjectured that standard quantum computation is more efficient than probabilistic computa- 
tion. The conjecture is supported by the ability to efficiently factor large numbers Q and simulate 
physics Q using quantum computers, by proofs that quantum computers are more powerful with 
respect to some black boxes 0, and by results showing exponential improvements in communi- 
cation complexity ||4|. 

To delineate the conjecture one can consider models of computation where the basic operations 
are multiplication of linear operators in a given set G. Each operator in G is associated with a 
complexity (e.g. the length of its name), so that the complexity of a product gig2 ... is the sum 
of the complexities of the gi. One can then ask questions about the complexity of calculating 
quantities like: 1. Computing the entries in a standard basis of a product. 2. Computing the trace 
of a product. When G is the a set of elementary quantum gates, the power of quantum computers 
is equivalent to being able to efficiently sample from a probability distribution with expectation an 
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entry of a product and variance 0(1) (see Q). The power of one-bit quantum computers Q is 
equivalent to sampling from a probability distribution with expectation the trace of a product and 
variance 0(2"'), where n is the number of qubits. 

A special case is when the set G is the group of operators normalizing the group generated by 
the Pauli matrices (bit flip, sign flip). For n qubits, this group has order 2'^^"^^ and plays a crucial 
role in encoding and decoding stabilizer codes |l6| and in fault tolerant quantum computation [7|. 
In lIHl it is shown that even when this group is extended by projections onto the logical states of 
qubits, the complexities of the two questions above are polynomial. Two similarly defined groups 
consist of the linear optics operators for fermions and for bosons. In both cases, the groups are Lie 
groups of polynomial dimension in the number of modes. (Modes play the same role as qubits in 
these systems). A few simulatability results were known for these groups. For example, for bosons, 
the orbit of the vacuum state under the linear optics operators consists of Gaussian states, for 
which many relevant quantities can be efficiently computed. Similarly, particle preserving linear 
optics operators applied to exactly one boson lead only to states that are equivalent to classical 
waves llQlfTOll. 

Recently, Valiant ifTDl demonstrated a set of products of operators (those definable by a class 
of "matchcircuits") for which the complexities of the first question and many of its generalizations 
are polynomial. Terhal and DiVincenzo lfT2ll realized that this set includes the unitary linear optics 
operators for fermions and that as a consequence, it is unlikely that it is possible to realize quantum 
computation in fermions by means of linear optics operators and particle detectors with feedback. 
They give a direct and efficient simulation of these operators based on fermionic principles. This 
result is at first surprising: In llT3l it was shown that with bosons, linear optics operators and particle 
detectors with feedback are sufficient for realizing quantum computation. The difference between 
fermions and bosons is explained by realizing that the effects of particle detectors are expressible 
as limits of non-unitary linear optics operators in fermions but not in bosons. As a result, the states 
achievable with fermionic linear optics operators and particle measurements are in the closure of a 
"simple" set. 

Since matchgate operators are non-unitary, one can ask what additional power is provided 
by Valiant's simulation of matchgates. Here it is shown that the two-qubit matchgates densely 
generate the monoid given by the closure of a group of extended fermionic linear optics operators 
in the Jordan-Wigner representation fT4l . This group defines the non-deterministic computations 
that can be physically realized with unitary linear optics operators and particle measurements. The 
equivalence of two-qubit matchgates and fermionic linear optics two-qubit operators generalizes 
to the set of simulatable matchcircuits introduced by Valiant. 
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2 Fermionic Linear Optics 

Let /, X^''\ Y^''\ Z^^'^ denote the identity and the Pauli operators acting on qubit k. Define Uk = 
^(1) ^(fc-i)f^(fc) (jj^ ^ ^Qj. If = x,Y. Then the Uk define a representation of fermionic 
mode operators. In particular, {X^ + iYk)/2 and {Xk — iYi)/2 represent the annihilation and 
creation operators for mode k. Let £i be the linear span of the identity together with the Uk for 
1 < k < n, where n is the number of modes (or qubits). The set Qi oi fermionic linear optics 
operators is the set of invertible matrices that preserve Ci by conjugation. That is, g E Qi iff for all 
A E Ci, gAg~^ E Ci. The terminology refers to the property that conjugation of an annihilation or 
a creation operator results in a linear combination of such operators. Let C2 be the set of products 
of two operators in Ci, so that C2 = C.\C\- The group Q2 of extended linear optics operators is the 
set of invertible matrices that preserve £2- Note that Q\EQ2. (In bosons, the analogous definitions 
lead to identical groups.) The group Q2 is considered to be "unphysical" for fermions, due to the 
presence of odd products of annihilation and creation operators. In Sect. |6lit is shown that Q2 is 
naturally viewed as a subgroup of Q\ for one more mode. 

The space £2 is a (complex) Lie algebra. It is spanned by the Pauli operator products given 
by /, f/fc, Z^^^ , and U^'')z^''+^^ . . . ^ith U,V E {X, Y}. The dimension of £2 is 

+ n + l. By considering general sums of Pauli products, one can check that if for every A E C2, 
[X, A] E C2, then X E £2- It follows that £2 is the Lie algebra of Q2. AH strictly quadratic (in £1) 
terms of C2, together with the identity also form a Lie algebra £3 of dimension — n, which is 
the Lie algebra of Qi. Physically, realizable operators are continuously generated from the identity. 
As a result, for the remainder of the paper, Qi is assumed to be given by the exponentials of £jQ 

In using (extended) linear optics operators for computation, one starts with the vacuum state 
|v.„) = |0 . . . 0)^ ,^ and applies operators in Qi (Q2) and measurements in the number basis |0)j^, 
for a mode. The outcomes of measurements are given by applying the measurement projections 
|0)fef=(0| = |(-^ + Z^''^) and 11)^^(11 = ^{1 — Z'^^'>). For standard computation, which projection 
"happens" is determined by the square amplitude of the result of applying it. For non-deterministic 
computation we can "choose" the outcome. In either case, analysis of the capabilities requires 
studying products of operators in Qi and the measurement projections. Let Si be the monoid given 
by the topological closure of Qi. 

If Q2 and measurements could be used for efficient faithful quantum computation, then the set 
of states Sn obtained with such operators from the n-mode vacuum state has to contain sufficiently 
large subspaces. That is, the 2'" dimensional state space of m qubits must be contained in Sn with 
n = 0(poly(m)). The following theorem makes this unlikely. 

'Without a proof that this assumption holds, it is possible that the groups studied here are only the component of 
the identity of the originally defined groups. 
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Theorem 1 The monoid generated by measurement projections and Q2 is contained in S2. 



Proof. This is a consequence of the fact that the measurement projections are limits of elements 
of G2- 

= lime-*^''Ve* (1) 

2 t^oo 



Since Q2 is a 2n^ + n + 1 -dimensional Lie group, Thm. [T] implies that iS2|v,i) is the closure 
of a small dimensional space. This suggests that S2 is not sufficiently strong for quantum compu- 
tation. The fact that the normalizer of the Pauli group together with standard measurements are 
insufficient ||8l follows in a similar way. That is, applying normalizer operators and projections 
onto stabilizer codes to the standard initial state always results in stabilizer states. 

Note that a similar result cannot be shown for bosonic linear optics operators with particle 
measurements. Only the projection operator onto the boson state of a mode is expressible as a 
limit of (non-unitary) linear optics operators. This provides an explanation of why efficient linear 
optics quantum computation is possible lfT3l . 

3 Matchgates and Linear Optics Operators 

In ifTTl , Valiant introduced a family of linear operators (called matchgates) acting on qubits. Match- 
gates are based on a graph theoretic construction. Valiant showed that under certain conditions, the 
coefficients of matrices defined by products of matchgates could be efficiently calculated. Match- 
gates acting on two qubits were shown to satisfy a set of 5 equations, the matchgate identities. If 
B is the matrix defined by a matchgate acting on two qubits, then the following are 0: 

Ml = (00|fi|00)(ll|fi|ll) - (10|fi|10)(01|fi|01) - (00|5|ll)(ll|fi|00) + (10|5|01)(01|5|10) 

(2) 

M2 = (10|5|00)(11|5|11) - (10|S|10)(11|S|01) - (11|5|00)(10|E|11) + (10|B|01)(11|5|10) 

(3) 

M3 = (01|5|00)(11|5|11) + (01|fi|01)(ll|fi|10) - (ll|fi|00)(01|5|ll) - (01|B|10)(11|5|01) 

(4) 

M4 = (00|fi|01)(ll|fi|ll) + (01|5|01)(10|fi|ll) - (00|fi|ll)(ll|B|01) - (10|5|01)(01|B|11) 
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(5) 

Ms = (00|5|10)(11|5|11) - (10|5|10)(01|5|11) - (00|5|11)(11|5|10) + (01|5|10)(10|5|11) 

(6) 

Let A^2 be the set of matrices B satisfying the identities Mi = and either (ll|i?|ll) 7^ or i? is 
diagonal. Valiant showed that these matrices are realizable by matchgates. 

Theorem 2 The closure of M.2 is S2for two modes (or qubits). 

Proof. The Lie algebra which densely generates ^2 is spanned by the 11 operators 

L = {II, XI, YI, ZI, ZX, ZY, XX, XY, YX, YY, IZ} (7) 

Here UV abbreviates f/^^V^^). One can check that for A G L \ {//}, A{YX) + {YX)A^ = 0: 
It suffices to note that if = A, then A anticommutes with YX, and if A'^ = —A, which 
is the case if A contains an odd number of F's, then A commutes with YX. (This property 
generalizes for arbitrary number of qubits, using the operator YXYX . . . instead of YX.) The 
identity A{YX) + {YX)A^ = can be re-written in the form {A(g)I + 1 (g) A)T = 0, where T is 
the antisymmetric vector 

T = |00)|11) - |11)|00) + |01)|10) - |10)|01). (8) 

This means that T is an eigenvector of the Lie group C generated byL©L = {Aigl + 1 ® A : A E 
L}. Note that C = {B ® B : B E ^2}- preserves antisymmetric vectors, so the statement that 
CT oc T is equivalent to R^CT = for all R antisymmetric such that R^T = 0. The dimension 
of such R is 5, and here is a basis: 



Ri = 


|00)|11) - 


|11)|00) - |01)|10) + |10)|01) 


(9) 


R2 = 


|00)|01) - 


|01)|00) 


(10) 


R3 = 


|00)|10) - 


|10)|00) 


(11) 


i?4 = 


|01)|11)- 


|11)|01) 


(12) 


i?5 = 


|10)|11)- 


|11)|10) 


(13) 



Define the expressions 

Ei = RjBT (14) 
Ef = T^BRi (15) 
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Since for two qubits C2 = C2, members B of Q2 satisfy the identities Ei = 0, Ef = 0. Because 
these identities are all derived from an eigenvector condition, the set of matrices B satisfying them 
is a closed monoid Q2 containing ^2- 
Using the equivalence 

{\ab)\cd)fB^B{\ef)\gh)) = {ab\B\ef){cd\B\gh) , (16) 



one can check that the following hold 



^1 + E^ 


= 4Mi 




(17) 


E4 


= 2M3 




(18) 




= 2M2 




(19) 


El 


= 2M4 




(20) 


El 


= 2M5 




(21) 


{ll\B\ll){Ei- El) 


= 4((01|5|11)M2 - 


(10|5|11)M3 + 


(ll|5|10)Af4 - (ll|5|01)Af5) 








(22) 


{ll\B\ll)E2 


= 2((01|fi|ll)A'/i - 


(00|fi|ll)Af3 + 


(01|fi|10)A'/4 - {01\B\01)M5) 








(23) 


(11|5|11)E3 


= 2((10|S|ll)Mi - 


(00|S|11)M2 - 


(10|E|01)A^5 + (10|S|10)A^4) 








(24) 


{ll\B\ll)El 


= 2((ll|5|01)Mi - 


(11|5|00)M4 + (10|E|01)M3 - (01|S|01)M2) 








(25) 


{n\B\n)El 


= 2((ll|fi|10)Afi - 


(ll|fi|00)Af5 - 


(01|5|10)Af2 + (10|fi|10)A//3). 



(26) 



Mathematica instructions to check the above relationships are included verbatim in Appendix lAl 

Since diagonal matrices trivially satisfy Ei = 0, Ef = (i > 1) and Ei — Ef = 0, the 
identities imply that M2 C Q'^. Let M'2 = {B e M2 : 7^ 0}. By directly solving for 

the entries of B other than in the first summand of the Mi, one can see that 7Vl2 ^^i 

analytically coordinatizable 11 complex dimensional manifold. The diagonal members of A^2 are 
in the closure of A^g- 

The identities also imply that the elements of Q2, and therefore those of S2, satisfy Mi = 0. It 
follows that the B e S2 with B diagonal or(ll|i?|ll) 7^0arein7Vl2- 

For invertible B, the identities = imply that B{XY)B^ = XXY for A 7^ 0. It follows 
that the tangent space at B is exactly that of Q2 at B. Consequently, and Q2 contain the same 
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invertible matrices satisfying (ll|i?|ll) ^ 0. It remains to show that these matrices are dense 
in both sets. For M.'^ it suffices to observe that for fixed ^ 0, there is an invertible 

i? G A^2' which implies that the determinant function is not null on this linearly defined subset. 
Hence the complement of the determinant's null set is dense. For the density property follows 
from the fact that the subgroup generated by XI and XX acts transitively on the basis states. 



4 Simulatable Matchcircuits 

Valiant showed that any composition of operators consisting of two qubit matchgates on the first 
two qubits and gates of the form e*(^*'°^^'''+^0 and e*(^''''^'''+^') is efficiently simulatable in the 
following sense: If i? is a product of m such gates, then many sums of squares or square norms of 
entries of B can be computed efficiently in m and n (the number of qubits). Let be the set of 
all products of the gates mentioned. 

Theorem 3 The closure of M. is S2. 

Proof. By definition and by Thm. El M. ^82- It suffices to show that the invertible operators 
in generate ^2- This can be checked directly by using the Bloch sphere rules for conjugating 
products of Pauli matrices by 90° rotations (e~*^'^/^) around other products ifTSll . For example, 
is obtained by conjugating Z^^^Y^"^^ with a rotation around X^'^^X^^\ The operator 
Z(^) is obtained by conjugating Z^^^ Z^'^^X^^^ with a rotation around Z^^^ Z^^^Y^''^^ The latter op- 
erator can be deduced similarly to the way Z^^^ Z^'^^X^^^ was obtained. Induction can be used to 
extend to arbitrarily many qubits. ■ 



5 Non-deterministic Computations 

A non-deterministic computation with fermionic linear optics consists of a sequence of linear op- 
erators and measurements, where one post-conditions on the measurement outcome in the sense 
that one multiplies the state by the appropriate projection operator. The outcome is not normalized. 
Let U be the implemented operator. The minimal quantities one wishes to compute efficiently are 
tr((v„|f/^(/ ± Z^''^)U\vn)), which give the relative probabilities of the outcome of a measurement 
on the /c'th mode. Suppose that implementable operators form a monoid and include the standard 
measurement projections. Since |v„)(v„| = Ilfcll-^ + and {I + Z) = {I + ZY{I + Z), it 
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suffices to be able to compute, for each implementable U, tr{WU) = J2ki \ Uki\'^- This motivates 
a definition that works for any monoid generated by elementary operators: An efficient simulation 
is defined to be an efficient algorithm for computing ix{WU) for an explicitly implemented (as 
a product of elementary operators) U . Efficiency is defined in terms of the implementation com- 
plexity of U . With this definition, Valiant demonstrated an efficient simulation of matchcircuits 
composed of certain matchgates. The purpose of this section is to discuss how that leads to an 
efficient simulation of a dense subset of the monoid S2 with naturally defined generators. 

An elementary fermionic gate is an operator of the form ae**^ = Q + with U one of the 
products of Pauli operators in £2 other than the identity. The coefficients q and r are required to be 
complex rationals with q ^ 0. Let d be the number of digits needed to denote these rationals. The 
description length of g + rU is il{2 log(n) + d), where the summand 2n is the description length 
of U, one of O(n^) many possible Pauli products. The elementary projection is the operator (/ + 
^(n+i))/2. It is implementable non-deterministically by post- selection on a particle measurement. 

The elementary fermionic gates can be realized in terms of the operators allowed in simulatable 
matchcircuits: Simply conjugate one of these operators by the appropriate sequence of 90° allowed 
operators. Note that the 90° operators are elementary if scaled by ^/2. The standard measurement 
projections are allowed in matchcircuits. It is therefore possible to take a product of elementary 
fermionic gates and projections, and efficiently express them using allowed matchgates. It follows 
that Valiant's algorithm can be used to efficiently simulate the monoid S2 generated by elementary 
fermionic gates and projections. The goal is to show that these operators densely generate ^2. 

Theorem 4 Except for a scale factor, the operators 0/^2 on n modes are implementable by first 
adjoining a mode in state |0),^^^, applying a sequence of unitary operators ofQ2for n + 1 modes 
and elementary projections and finally discarding mode n + 1. 

Proof. It suffices to show that e*^'"' with real t is implementable up to a scale. This follows from 
the observation that other real exponentials of Pauli operators are conjugates of e*^*"' by unitary 
operators, and these together with unitary operators generate ^2- 
To implement e^*"*^'"' realize the following sequence of operators: 

1. Adjoin |0),^_|_^ (if that hasn't already been done). 

2. Apply e**(^*"'^*"^''+'*^*"'^'"^'')/2 

3. Project mode n + 1 with (/ + Z^'^^^^)/2, which returns mode n + 1 to its initial state, or 
results in 0. 
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To see how this works, apply it to a|0)„ + Step 2 is a partial swap with a phase and 

results in a|0)„ + /3{cos{s)\l)j0)^_^-^ + z sin(s)|0),Jl),^_|_J. The elementary projection results in 
(a|0)^ + cos(s)/3|l)„)|0)„_|_j. It follows that the effect is the same as applying a scalar multiple of 
g-in(cos(s))z(")/2 rpj^g other sign in the exponent can be obtained by replacing step 2. with: 

2' Apply e««(^*"'-^^"+^'-'*^'"'^*"+^')/2 



Corollary 5 The closure of £2 is ^2 ® (/ + Z("+i))/2. 

Proof. This follows from Thm. |4| and the fact that the elementary rotations e**^ for U a Pauli 
product densely generate all such rotations. (See, for example, lfT6ll .) ■ 

It can be seen that the ability to efficiently simulate non-deterministic computation as defined 
above leads to an efficient simulation of a quantum computation with measurements and future op- 
erators conditioned on the measurement outcomes. The method is described in lfT2l and basically 
consists of simulating, at each step, the random measurement outcome, using a calculation of the 
conditional probability distribution. 

A potentially easier problem then efficient simulation of a monoid is to determine, for an im- 
plemented U, whether U = 0. Observe that if it was possible to use ^2 with elementary generators 
to efficiently and faithfully realize quantum computation, then the zero-test algorithm can be used 
to efficiently solve problems in polynomial quantum non-deterministic time as defined in ifTTll . 
In ifTSi it was shown that this is hard for the polynomial hierarchy. 

6 Identifying the Lie Algebras: Qi is General 

Let £2 be the set of trace zero members of C2. The adjoint action of Q2 on C2 permits representing 
members of Q2 as (2n^ + n) x (2n^ + n) matrices. The representation is faithful up to scalar multi- 
ples, because C2 algebraically generates all operators on the n qubits. This means that products of 
elementary operators can be efficiently computed in the representation. The reverse procedure, i.e. 
finding a decomposition of a represented operator in terms of a product of exponentials of Pauli 
products is also possible, though less obviously so. For this purpose it is more useful to recognize 
C2 as the Lie algebra S02n+iC and work in the fundamental representation. One way to recognize 
C2 is to realize that it is (isomorphic to ) a subalgebra of C'2 for one more qubit. The mapping is 
accomplished by modifying the members of the form Uk by multiplying with X^^\ This makes 
the operators strictly quadratic for fermionic modes 0, . . . , n (in this order). Then observe that 



9 



adjoint action on Ci is the fundamental representation of 502n+2C- The algebra can now be 
identified. Incidentally, this construction shows that in a sense C2 is no more general then C'2 de- 
spite appearances. This together with the results of the previous section implies that the simulation 
algorithm of Terhal and DiVincenzo ifTlll can be used to simulate ^2 with the same generality as 
Valiant's. 

Here is the direct way to identify £3 ^ Lie algebra: In the fundamental representation 
S02n+iC is spanned by the antisymmetric matrices Sij = — for < i < j < n. 

The identification is made via the correspondences 





^Xk/2 - 




(27) 




iYk/2 - 


So{n+k) 


(28) 




iZ(^V2 - 


Sk{n+k) 


(29) 




. .X(^V2 - 


S{n+l)k 


(30) 




. .r(^V2 - 


S{n+l)(n+k) 


(31) 


,y{0^(/+i) . 


_X(^-)/2 - 


-Slk 


(32) 


,y{0^(/+i) . 


. - 


—Sl{n+k) 


(33) 



This identification of C2 permits efficiently representing a product of elementary operators as 
a {2n + 1) X (2n + 1) matrix. Let A be a matrix thus obtained. Then A^A = I, and this identity 
characterizes the Lie group generated by the Ski- The process of representing a matrix satisfying 
A^A = / as a product of elementary operators is straightforward by using a variant of Gaussian 
elimination to represent A as a product of O(n^) matrices of the form e**"*'-' = (J+s^,;) — cos(t)s^; + 
i sm{t)ski (t may be complex). By using conjugation rules by 90° rotations, one can then expand 
this into a 0{n^) product consisting only of operators that are allowed for Valiant's simulatable 
matchcircuits. 

7 Concluding Comments 

It is true that bosons can be represented by paired fermions. So why does this not lead to an ef- 
ficient realization of quantum computers by using this representation together with techniques for 
bosonic linear optics? One answer is that the bosonic linear optics operators in this representa- 
tion correspond to Hamiltonians that are quartic in the annihilation and creation operators and are 
therefore not in £2- It is in fact not hard to see that adding to £2 only the Hamiltonian Z^^^Z^^^, 
the Lie algebra generated contains all products of Pauli matrices and so generates all invertible 
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matrices lfT9l . 



Suggested problem areas for future investigations: 

1 . Determine the complexity of efficiently simulating representations of the three families of 
simple complex Lie groups. Is the complexity polynomial in the dimension of the groups? 

Notes: 

The results of Valiant, Terhal and Divincenzo and this paper show that the answer is 
"yes" for one family of representations. 

The answer might depend on the choice of generators and elementary operators. The 
fundamental representation of each such group can be used to make a reasonably natu- 
ral definition. 

Which projectors in a representations are to be assumed as elementary operators? They 
should be in the closure of the Lie group. 

Semisimple Lie algebras can be analyzed in terms of their simple parts. What about 
non-semisimple ones? 

2. What finite monoids of operators are efficiently simulatable? 
Notes: 

Again, the choice of generators may be crucial, and it is desirable that it is "natural" in 
some sense. 

The monoids associated with n-ary stabilizer codes via the appropriate normalizer are 
efficiently simulatable in terms of the number of systems used. 

Is the stabilizer code example naturally generalizable? 

3. Problem areas 1. and 2., but for efficiently determining whether a product of generators is 
zero. Is this sometimes strictly easier to do? 

4. Find a group or monoid of operators where the probabilistic behavior of a (quantum) com- 
putation is efficiently simulatable, but the non-deterministic behavior is not. 

Notes: 

It is necessary to define what is meant by "probabilistic" behavior. The one case where 
an interpretation is readily available is if the group is unitary and the initial state as 
well as standard measurements are provided. For a monoid, one approach is to allow 
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as measurements some or all partitions of unity definable by its operators. The monoid 
should be (densely) generated by its unitary operators and projections associated with 
measurements. 
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A Checking the Matchgate Identities 

(* Mathematica notes. *) 

(* Useful rules: *) 
Unprotect [Dot] ; 

Dot [tensor [ a_, b_] , tensor [ c_, d_] ] = (a . c) * (b . d) ; 
Dot[-a_,b_] = -(a.b); 
Dot[a_, -b_] = -(a.b); 
Dot[-a_, -b_] = (a.b); 

(* For obtaining the equation for the transpose: *) 
trnsprls = {b[c_].k[d_] -> b[d].k[c]}; 

(* For obtaining the equation for the conjugate by XX: *) 
xxrls = {xOO->xll, x01->xlO, xlO->x01, xll->xOO}; 
(* Swapping: *) 

swprls = {x01->xlO, xlO->x01 } ; 
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Iswprls = {b [xOl] ->b [xlO] ,b [xlO] ->b [xOl] } ; 

(* Conventions: 

* b[xab] stands for $\bra{ab}$, k[xab] for $\ket{ab}$. 

* Quadradic expressions for a matrix B are expressed 

* $\trace(X (B\tensor B) ) $ with X in the appropriate 

* tensor product space. X is given for various expressions. 

* This way the expression (b [xO ] . k [xO ] ) * (b [xll ] . k [xO 1 ] ) 

* refers to the product $ \bra { } B\ket { } \bra { 11 } B\ket { 1 } $ . 
*) 



(* Matchgate expressions: *) 



Ml 








[xOn 1 k 


r xi 1 1 1 + 

2Si _L _L J J ^ 






— "h f^i Ti Q o "K" 




1 "h ^ Ti Q O "K" 


rk r vi n 1 


k r vDl 1 1 

jS- L ^ U -L J J 






— "h ^ Ti Q O "P" 




1 "h Ti Q O "K" 


rk rvi 1 1 


k r vDD 1 1 

jS- L ^ u u J J 










1 -h ^ f-i o o 7" 


rk rxDi 1 


k r xi n 1 1 


r 




— j~ ^ T\ C! (~\ Y' I 






r vDD 1 k 


r vi 1 1 1 + 

2i J. J. J J T 








rhrxini hrxiii 


1 -h pi f-i o o 7" 


rk rxi n 1 


k rxDI 1 1 


+ 




- tensor 


.bLxllJ ,bLxl(JJ 


] . tensor 


.k [xUU J , 


k [ xi 1 J J 


+ 




+ tensor 


[b [xlO] ,b [xll] 


] . tensor 


[k [xOl] , 


k [xlO] ] 


r 


M3 


= tensor [b 


[xOl] ,b [xll] ] . 


tensor [ k 


[xOO] , k 


[xll]] + 






+ tensor 


[b [xOl] ,b [xll] 


] . tensor 


[k [xOl] , 


k [xlO] ] 


+ 




- tensor 


[b [xll] ,b [xOl] 


] . tensor 


[k [xOO] , 


k [xll] ] 


+ 




- tensor 


[b [xOl] ,b [xll] 


] . tensor 


[k [xlO] , 


k [xOl] ] 


r 


M4 


= tensor [b 


[xOO] ,b [xll] ] . 


tensor [ k 


[xOl] , k 


[xll]] + 






+ tensor 


[b [xOl] ,b [xlO] 


] . tensor 


[k [xOl] , 


k [xll] ] 


+ 




- tensor 


[b [xOO] ,b [xll] 


] . tensor 


[k [xll] , 


k [xOl] ] 


+ 




- tensor 


[b [xlO] ,b [xOl] 


] . tensor 


[k [xOl] , 


k [xll] ] 


r 


M5 


= tensor [b 


[xOO] ,b [xll] ] . 


tensor [ k 


[xlO] , k 


[xll]] + 






- tensor 


[b [xlO] ,b [xOl] 


] . tensor 


[k [xlO] , 


k [xll] ] 


+ 




- tensor 


[b [xOO] ,b [xll] 


] . tensor 


[k [xll] , 


k [xlO] ] 


+ 




+ tensor 


[b [xOl] ,b [xlO] 


] . tensor 


[k [xlO] , 


k [xll] ] 


r 



(* Check: 

M3 - (M4/ . trnsprls) 
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* = Q * 

M2 - (M5/ . trnsprls) 

* = Q * 



(* Lie expressions: *) 
T = tensor [k [xOO] , k [xll] ] 
tensor [k [xOl] , k [xlO] ] 
Rl = tensor [b [xOO] ,b [xll] 
tensor [b [xlO] ,b [xOl] 
R2 = tensor [b [xOO] ,b [xOl] 
R3 = tensor [b [xOO] ,b [xlO] 
R4 = tensor [b [xOl] ,b [xll] 
R5 = tensor [b [xlO] ,b [xll] 



tensor [k [xll] , k [xOO] ] + 
tensor [k [xlO] , k [xOl] ] ; 

- tensor [b [xll] ,b [xOO] 

- tensor [b [xOl] ,b [xlO] 

- tensor [b [xOl] ,b [xOO] 

- tensor [b [xlO] ,b [xOO] 

- tensor [b [xll] ,b [xOl] 

- tensor [b [xll] ,b [xlO] 



El = Distribute [Rl . T] ; 
ETl = El/ .trnsprls; 
E2 = Distribute [R2 . T] ; 
ET2 = E2/ .trnsprls; 
E3 = Distribute [R3 . T] ; 
ET3 = E3/ .trnsprls; 
E4 = Distribute [R4 . T] ; 
ET4 = E4/ .trnsprls; 
E5 = Distribute [R5 . T] ; 
ET5 = E5/ .trnsprls; 
(* Clneck: 

Simplify [El+ETl - 4*M1] 

* = Q * 

Simplify [E4 - 2*M3] 

* = Q * 

Simplify [E5 - 2*M2] 

* = Q * 
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Simplify[ (b [xll ] . k [xll ] ) * E2 - 
2* ( 

(b [xOl] .k [xll] ) *M1 + 

- (b [xOO] .k [xll] ) *M3 + 
(b [xOl] .k [xlO] ) *M4 + 

- (b [xOl] .k [xOl] ) *M5 



Simplify[ (b [xll ] . k [xll ] ) * E3 - 
2* ( 

- (b [xOO] .k [xll] ) *M2 + 
- (b [xlO] .k [xOl] ) *M5 + 

(b [xlO] .k [xlO] ) *M4 + 

(b [xlO] .k [xll] ) *M1 



Simplify [ (Ml/ . trnsprls) - Ml] 

* = * 
* 

Simplify [ (M2/ . Iswprls) - M3] 

* = Q * 

Simplify [ (E2/ . Iswprls) -E3] 

* = * 

Simplify [ (b [xll] .k [xll] )* (El - (El /. trnsprls ) ) - 
4* ( 

b [xOl] .k [xll] *M2 + 
-b [xlO] .k [xll] *M3 + 
-b [xll] .k [xOl] *M5 + 

b [xll] .k [xlO] *M4 
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= 0* 

This confirms the identites claimed in the text. 
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